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Abstract. A positive experimental measurement of the nuclear Schiff moment would 
have important implications for physics beyond the standard model. To aid in the 
interpretation of a proposed experiment to measure the nuclear Schiff moment of ^°''Pb 
in the ferroelectric PbTiOa, three-dimensional Hartree-Fock calculations have been 
performed to model the local electronic structure in the vicinity of the Pb nucleus. 
The energy shift due to the Schiff moment is found to be a factor of 2 smaller in 
comparison to existing estimates. 
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1. Introduction 

The search for the permanent electric-dipole moment (EDM) of quantum particles has 
been of continued interest for over 40 years since the discovery of the violation of the 
combined symmetry of charge conjugation (C) and parity (P) in the decay of the 
meson [1]. By the CPT theorem, the existence of CP- violation also implies that time- 
reversal (T) symmetry is violated [2]. P- violation, together with T-reversal asymmetry 
gives rise to a permanent electric-dipole moment of a quantum system in a stationary 
state. Therefore ongoing searches for EDM's of elementary particles, nuclei, atoms and 
molecules are important for studies of fundamental symmetries [3] and they provide 
important constraints on theories that attempt to go beyond the standard model 

Currently, no EDM has been experimentally observed. Experiments on 
paramagnetic atoms provide the best upper limit on the electron EDM, with the most 
stringent limitation coming from experiments with atomic thallium [5] and the YbF 
molecule [6] 

Tl: 4 < 1.6 X lO^^^ecm. 

YbF: 4 < 1.05 X 10~^^ecm . (1) 

For diamagnetic atoms, the most EDM sensitive experiment has been performed with 
mercury vapour, ^^^Hg, giving an upper limit for the EDM of the ^^^Hg atom of [7], 

rf(^^^Hg) < 3.1 X lO-^^ecm. (2) 

For diamagnetic atoms the major contribution to the electron dipole moment is from 
the nuclear Schiff moment (NSM), S. The NSM can be defined by the P- and T-odd 
electrostatic potential [8], 

(/p(r) = 47r(S ■ V)5(r). (3) 

The result of [7] has been interpreted with the aid of calculations P HO] to yield an 
upper limit of the Schiff moment of the ^^^Hg nucleus of, 

^(^^^Hg) < 0.8 X 10-'^ea|, (4) 

where is the Bohr radius. 

Currently, there is intense interest in exploring physical systems that potentially 
can provide orders of magnitude increases in sensitivity to CP violating effects. As 
first discussed by Shapiro [H], condensed matter systems are particularly promising in 
this regard. The electron's EDM is aligned with the spin of the electron, and therefore 
its magnetic moment. It follows that in a compound with uncompensated spins, the 
application of an external electric field will align the EDM's and thus align the magnetic 
moments, giving rise to a macroscopic magnetization. Experiments then aim to measure 
the reversal of the magnetization by reversing the direction of the external electric field. 
Recent efforts have focused on gadohnium garnets [121 ESI [El CS] and more recently on 
the ceramic Euo.sBao.sTiOs [161 [HI HH [19] . 

Similarly to the atomic case, condensed matter systems without unpaired electrons 
should also be sensitive to the NSM. Ferroelectric PbTiOs was proposed for searches 
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of time reversal symmetry violations enhanced by the internal electric field [2U] . 
Experiments with this compound can potentially lead to improvement of sensitivity 
to the sNSM by several orders of magnitude in comparison with existing results 
[2111221 [231 [Ml 125]. 

The NSM can be detected experimentally either by a measurement of the 
macroscopic magnetization induced by an electric field or by nuclear magnetic resonance. 
The lead nuclei will interact with the large internal electric field in the ferroelectric 
leading to a P, T-odd energy shift. Here, we perform three-dimensional Hartree-Fock 
calculations in order to model the local electronic structure of PbTiOa around the Pb 
nucleus and determine the energy shift caused by the NSM. 

2. Model calculations 

PbTiOs is an ionic crystal made up of Pb^+, Ti^+ and ions. In the ferroelectric 
phase, PbTiOg is tetragonal with c/a = 1.065, a = 3.902A, c = 4.156A [26], with Fh'^+ 
and Ti^+ displacements of 0.47A and 0.30A respectively [27] . The ion displacements 
are in accordance with the Schiff theorem [28]. The lattice relaxes so as to screen the 
strong internal ferroelectric field such that the average field that acts on each charged 
particle in the compound is zero. However, the finite size of the ions plays an important 
role, with oxygen electrons (s and p-wave) penetrating into the core of the Pb ion in 
the PbOi2 cluster (Fig. [1]). Due to the displacement of the Pb ion, this creates an 
electric field gradient at the Pb nucleus. The NSM interacts with this gradient, leading 
to an energy shift. The energy shift can be calculated using second order perturbation 




Figure 1. The Pb^+ ion surrounded by the first coordination sphere of 12 ions. 
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theory. In Ref. [13], using a semiclassical approximation for the electron wavef unctions, 
a formula for the energy correction due to the NSM was obtained, 

ME, ^ {\R.n + ^^3/2) (5) 



PI 



)3/2 ' 3-^^/^; ^^e^i 



where Ae is the energy shift caused by the NSM S [29], X is the ion displacement, 
Eq = 27.2 eV, e = |e| is the charge of the electron, Z = 82 is the nuclear charge 
of Pb, Us and Up are effective principal quantum numbers and b = 16/3/a/3. /3 is a 
parameter describing lead orbitals partially occupied by oxygen electrons. R1/2 and 
R3/2 are relativistic enhancement factors [3l [13] given by, 

A 271/2-2 

-Kl/2 ~ 



[r(27i/2 + 1)]^ 

P , 487V24-^-""^-^ f6l 

r(27i/2 + i)r(273/2 + i)' 

where 71/2 = y/T^^'Zh?, 73/2 = a/4 — Z'^a^, a is the fine structure constant, T{x) is 
the gamma function and Xq = {2ZrQ/aB), with rg the nuclear radius. 

An understanding of the order of magnitude of this energy shift can be gained by 
considering that the energy shift can be approximated by the product of the internal 
electric field in ferroelectric PbTiOs, Eint — 10^ V/cm, and the EDM of the Pb^"*" ion, 
d(Ph'^~^)/{ecm) ^ 10~^S'/(ea^), where we have assumed that the EDM of the Pb^"*" ion 
is similar to that of the Hg atom. This gives, Ae ^ E^diPh^^) ^ lO^^eV [22]. 

In Ref. [21], Eq. (|5]) was applied to calculate the energy shift using a value of 
/3 = —0.29 estimated previously for GdOg clusters [T3] . 

This gave an estimate for the energy shift of, 

S 

Ae ^ -1.1 X lO^^eV. (7) 



The method used to calculate /3 involved matching a linear combination of the 
electron wavefunctions of the 2po- electrons of possessing the correct cubic symmetry, 
with single particle 6s and 6p states of the central Gd^"*" ion that can be thought of as 
being occupied by the oxygen electrons. The matching is done at an intermediate 
distance in the cluster between the ions and the central Gd'^"'" ion. 

Aside from the desirability of using a value of /3 calculated specifically for the PbOi2 
cluster, the fact that this cluster is closed shell and free from complications arising from 
open shell effects allows a more accurate method of calculating /3 to be used in order to 
benchmark the earlier estimates. By relating /3 to a shift in the electron density around 
the Pb nucleus, the calculation of /3, and thus the energy shift due to the NSM, is 
reduced to determining the electron density of the PbOi2 cluster. This allows the well- 
developed theoretical apparatus of quantum chemistry to be applied. Once a ground 
state wavefunction is obtained via the Hartree-Fock approximation, the density shift 
around the Pb nucleus can be simply extracted. 

First we relate the density shift to /3. Ap{z) will be defined as the electron density 
shift around the lead nucleus due to the displacement of the lead ion from the center 
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of the PbOi2 cluster, z is the distance from the lead nucleus in the direction of the 
displacement. This can be represented by, 

Ap(.) = fe<Z£). (8) 

It is possible to make a simple estimate of the density shift as follows. Near the Pb^"^ 
ion the electrostatic potential for electrons is dominated by that of the central Pb^^ ion. 
Therefore the wavefunctions of electrons in this region can be well described by partially 
occupied Pb^"*" orbitals. The interaction of electrons with the NSM of Pb is highest near 
the Pb nucleus. The largest contribution to this interaction comes from s— and p-wave 
states that can penetrate close to the nucleus. First we let IS") and \P) be a combination 
of the Pa states of the PbOi2 cluster that have s— and p— wave symmetries with respect 
to the centre of the cluster. Then these states can be represented as partially occupied 
7s and 6p orbitals of the central lead ion. 

\S)^l3s\7s), |P,)^/3p|6p,), (9) 

where i = {x, y, z}. When the Pb ion is shifted from the centre of the cluster, the IS*) 
and \P) states no longer have exact s and p— wave symmetries with respect to the lead 
ion but are instead a mixture of the two. If we assume the linear approximation, then 
the admixture of states of different parity due to the displacement X of the Pb^"*" ion is 
of the form, 

1^) ^(3s\7s) + Y.(3's—m 

i 

|P.)^/3p|6p,,) + E/3;— (10) 

i «B 

The change in the density due to the displacement of the Pb^"*" ion along the z direction 
is then, 

Ap{z)=4/3—ij^,{z)i:,p{z), (11) 
as 

where (3 = {PsP's + PpP'p) and a factor of 2 comes from the double occupancy of each 
spin orbital. 

In order to obtain the electron density in the vicinity of the lead nucleus, we adopt 
a simple model of the electronic structure with only the Pb^"*" ion and the 12 ions 
in the first coordination sphere of the Pb^"*" ion being included ab-initio. The 12 
ions are placed at the mid-points of the edges of a cube of length 4A. The Pb^"*" ion 
is displaced by a distance X from the centre of the cube and 8 Ti*^"^ ions (see below) 
are represented by point charges displaced from the corners of the cube by 0.3X/0.5. 
Varying the Pb ion displacement will then also change the Ti ion displacement, this will 
model a change in the strength of the internal ferroelectric field. Point charges are also 
added to ensure that the system is neutral. These charges are placed at the positions 
of the 24 ions surrounding the PbOi2 cluster. 

The GAMESS-US quantum chemistry code [30] is used to perform Hartree-Fock 
calculations for PbOi2. The WTBS basis [3ll[32] set is used for all calculations presented 
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here. A diffuse sp shell is also added to the Pb^"*" and orbitals in order to better 
represent the ionic bonds present in this system. A trial wavefunction is built from 
Pb^"*" and Hartree-Fock wavefunctions, where in both cases the ions are enclosed 
by a cubical array of point charges. For this is necessary in order to obtain a stable 
ground state in the Hartree-Fock approximation. This trial wavefunction is then used 
as a starting approximation to obtain a converged Hartree-Fock wavefunction for the 
PbOi2 cluster. Using a converged Hartree-Fock wavefunction for the PbOi2 cluster, the 
density shift around the Pb nucleus can be calculated, and then via Eq. ( ITT]) used to 
obtain a value of /3. The energy shift due to the NSM is then found via Eq. ([5]). 

3. Results 

In order to test the various assumptions made in the present model, a number of 
calculations have been carried out. Firstly, the point charges representing the Ti ions 
are treated here as an adjustable parameter in order to assess the dependence of the final 
result on the crystal structure, with calculations carried out for a range of Ti charges. 
Secondly, the assumed linearity of the density shift is tested by carrying out calculations 
for a number of Pb displacements. 

First, to illustrate how well Eq. ( ITT]) performs in capturing the essential physics 
of the density shift, Ap as a function of the distance from the Pb nucleus is plotted in 
figure [2] for a Ti charge of -1-3.5. We match the numerical density shift to the analytical 
expression of Eq. (|TT|) using atomic Hartree-Fock wavefunctions for the 6p and 7s states 
of Pb^"*", treating /3 as an adjustable parameter. The matching is done at points up 
to the first node, and then an average value of /3 is taken. As is clear, the analytical 
expression agrees well with the numerical results. 




0.05 0.1 0.15 0.2 



Figure 2. /^p{z) for X=0.5A, q = +3.5. Solid line, GAMESS-US; dashed line, Eq. 
(fll|) with P = -0.125. 
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Next, the validity of the hnear approximation for the density shift is investigated. 
Figure [3] plots Ap at distances from the Pb nucleus corresponding to the first minimum 
and maximum of Ap, as a function of the Pb^"'' ion displacement. From the figure, it is 
seen that a Ti charge of 3.5 gives the best agreement with the linear approximation out 
to X = O.SA, with some deviation from linearity observed for charges of 3.0 and 4.0. 



0.1 
0' 

1-0.2 
-0.3 
-0.4 

"^•^0 0.1 0.2 0.3 0.4 0.5 0.6 

X (Angs) 

Figure 3. Density shift as a function of the Pb ion displacement X, (with Ti charges 
displaced by 0.3X/0.5) for varying Ti charges q. Dashed line, q = +3.0; sohd line, 
q = +3.5, dot-dashed, q ~ +4.0. Lower lines at z = 0.004A, upper lines at z = 0.02A 

Finally, /3 as a function of the Ti charge is plotted in figure |4]for a Pb^"*" displacement 
of X = 0.5A. It is seen that for Ti charges between +3 and +4, /3 ranges from -0.1 to 
-0.125, demonstrating that /3 does not depend strongly on the crystal structure. This 
compares to a previous estimate of /3 = —0.29 [211 122] used for PbTiOs that was orginally 
derived for the garnet structure of a Gd'^+ ion enclosed by eight ions [13]. Taking 
a value of /3 ~ —0.125 for a Ti charge of +3.5 as an upper limit for Eq. ([5]) is 
evaluated for a Pb displacement of 0.47A, with effective principal quantum numbers of 
the 6p and 7s states of Pb^^, derived from Hartree-Fock energies of 1.028 and 1.417 
respectively [33]. The nuclear radius, tq, of ^"^Pb is estimated by the approximation, 
ro = CA^^^, where A is the atomic mass number and C = 1.25 x 10~^^ m. This gives 
ro — 7.39 X 10~^^ m. The resulting energy shift is then found to be, 

Ae ~ -0.66 X 10^ ^"^" eV = -0.59 x 10*^^ eV. (12) 
aBea% ea% 

This is approximately two times lower than the existing estimate, see Eq. (|7|). 
4. Conclusions 

The work in this paper is part of a continuing theoretical and experimental programme 
in aid of experimental investigations of CP violating effects in condensed matter systems. 
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Figure 4. /3 as the Ti charge is varied for X ~ O.SA. 



A Hartree-Fock calculation of the PbOi2 cluster enabled the density shift around the Pb 
nucleus due to the penetration of the oxygen electrons into the Pb ion to be calculated. 
From this an energy shift was calculated that is almost a factor of 2 lower than a previous 
estimate. 

There are two areas in the current work that merit further investigation. Firstly, 
the effect of the potential due to the crystal field needs to be properly taken into account. 
A possible solution to this issue is via the CRYSTAL code [Hll [35] . This particular code 
uses linear combinations of atom centred Gaussian functions to perform both Hartree- 
Fock and a variety of density functional theory calculations. This would allow a good 
description of both the crystal structure and also the penetration of the O electrons into 
the Pb ion. An existing calculation using the CRYSTAL code has provided accurate 
structural and electronic properties of PbTiOs [36] . 

Secondly, relativistic effects due to the heavy Pb nucleus need to be taken into 
account in order to correctly describe the behaviour of O electrons close to the Pb 
nucleus. As a first step, this can be done by solving the Dirac-Fock equations for the 
PbOi2 cluster. A challenge here will be in converging the self-consistent 4-component 
Dirac-Fock equations for such a large system. The BERTHA [37] or DIRAC codes [3H] 
can be investigated for this purpose. 

These two phenomena that must be included in any complete treatment of the 
problem, namely the potential due to the crystal lattice and relativistic effects due to 
the high nuclear charge of the Pb nucleus, make this a highly demanding theoretical 
and computational problem. 

A positive experimental finding of the nuclear Schiff moment would be of great 
importance to the entire physics community. Theoretical input into the design and 
interpretation of future experiments will be crucial in this effort. As such, we hope that 
the present model calculations will prompt renewed theoretical interest in this field. 
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